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ABSTRACT 



43 ■ 

O ' Context. Determinations of beryllium abundance in stars, together with lithium, provide a key tool to investigate the so far poorly understood 
extra-mixing processes at work in stellar interiors. 



^2 Aims. We measured Be in three open clusters, complementing existing Be surveys, and aiming at gathering a more complete empirical scenario 
i i ' of the evolution of Be as a function of stellar age and temperature. 

Methods. We analyzed VLT/UVES spectra of members of NGC 2516, the Hyades, and M 67 to determine their Be and Li abundances. In the 
t-H ■ first two clusters we focused on stars cooler than 5400 K, while the M 67 sample includes stars warmer than 6150 K, as well as two subgiants 
^ ' and two blue stragglers. We also computed the evolution of Be for a 0.9 M G star based on standard evolutionary models. 

, Results. We find different behaviours for stars in different temperature bins and ages. Stars warmer than 6150 K show Be depletion and follow 
■ a Be vs. Li correlation, while Be is undepleted in stars in the ~ 6150-5600 K range. NGC 2516 members cooler than 5400 K have not depleted 
| any Be, while older Hyades of similar temperature show some depletion. Be is severely depleted in the subgiants and blue stragglers. 
. . Conclusions. The results for warm stars are in agreement with those of previous studies, supporting the hypothesis that mixing in this tempera- 
•/"") ' ture regime is driven by rotation. The same holds for the two subgiants that have evolved from the "Li gap". This mechanism is instead not the 
dominant one for solar-type stars. Be depletion of cool Hyades cannot simply be explained by the effect of increasing depth of the convective 
zone. Finally, the different Be content of the two blue stragglers suggests that they have formed by two different processes (i.e., collisions vs. 
binary merging). 
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1. Introduction (those including convection only). In particular, observations 

of Li in open clusters of different ages have shown that these 

Beryllium (Be) is easily burned in stellar interiors at tempera- . , , . .. , . . ,,,„> . .. t 

J v ' J r stars deplete Li during the main sequence (MS), in spite of 

tures above ~ 3.5 MK; as a consequence, any process that is t , . . .. , , , , , 

M J r the fact that their convective zones do not reach deep enough 

able to mix surface material with internal Be-depleted material . „ . . . . c T . , . , 

r in the stellar interior to meet the zone of Li destruction (e.g., 

results in surface Be depletion. Since the more fragile lithium n ,. , lonn/c i a t , u • •> tU t » t mnt- 

r b Randich 2006 and references therein); the factor of ~ 100 Li 

(Li) is destroyed at temperatures above ~ 2.5 MK, measure- , , .. , . „ .. , £ . , , , 

v ' J r depletion measured in the Sun itself indeed represents a long 

ments of both elements in the same star allow us to investi- . in* c „ c i n «. • t • • u 

standing puzzle. Moreover, a factor of ~ 10 scatter in Li is ob- 

gate how deep mixing extends in the stellar interior and to put , . , , . „. .. , . £n /c .. 

b r b r served in the solar-age, solar metalhcity cluster M 67 (Spite et 

constraints on the various extra-mixing processes introduced , i lr . ~ i ^ , T , . , i lr . oo l „ . . , , i lr . r .- l T 

fe F al. 1987 1 Garcia Lopez et al. 1988 ; Pasquini et al. 1997 ; Jones et 

in evolutionary models to explain the complex Li behavior. , i lr , rir .i . ,. . . ... . ,. . 

1 ±- r al. 1999), indicating that otherwise similar stars (i.e., stars with 

Several unexpected results related to the depletion of Li in stars , , , , ,,. .. s , a . , , 

r r the same age, temperature, and metalhcity) may be affected by 

in fact have been found, which reveal our limited understanding . f1 - , , .• • • . 

6 different amounts of Li depletion, again in disagreement with 

of the physics acting in the interiors of stars. . , , , , rr U . ,■ ■ ■ ,. . T • , , .■ 

r J b standard models. This dispersion implies that Li depletion must 

Focusing on late-F and early G-type stars, different empir- fee affected fey an additional parameter besides mass, age, and 

ical evidences contradict the predictions of standard models chemical composition. 

Send offprint requests to: S. Randich, email: randich® arcetri.astro.it As to warmer early-F stars, since the 80s' Boesgaard and 

* Based on observations collected at ESO-VLT, Paranal collaborators pointed out the presence of a so-called "Li gap" 

Observatory, Chile, Programme numbers 65.L-0427, 68.D-0491, or "Li dip" among hot MS stars in the Hyades cluster, namely 

69.D-0454 a sharp decrease of Li in a narrow (~ 300 K) effective tempera- 
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ture (T e ff) range around 6600 K (Boesgaard & Tripicco fl986l l. 
The gap is present in several other clusters of the Hyades age 
or older (Balachandran 1995 and references therein), is weak 
or absent in the 120 Myr Pleiades (Pilachowski et al. 119871 
Boesgaard et al. 119881 ), but has been detected in the slightly 
older M 35 (160 Myr -Steinhauer & Delivannis 120041 1. sug- 
gesting that it develops at ages below ~ 200 Myr. 

In order to explain the unexpected MS Li depletion of 
F and early-G type stars, a number of models and extra- 
mixing processes have been proposed, which include mass 
loss (Swenson & Faulkner [T9921 , diffusion (Michaud [19861 
Chaboyer et al. 119951) . meridional circulation (Charbonnel & 
Talon 119991 and references therein), angular momentum loss, 
and rotationally driven mixing (Schatzman & Baglin 1 199 1 1 
Pinsonneault et al. 119921 Deliyannis & Pinsonneault 1997), 
gravity waves (Garcia Lopez & Spruit 119911 Montalban & 
Schatzmann |2000), tachocline (Brun et al. 119991 Piau et al. 
2003), and combinations of these (Charbonnel & Talon 2005). 
Note that models including rotational mixing would in prin- 
ciple be able to explain the scatter il Li in M 67, with initial 
angular momentum being the additional parameter affecting Li 
depletion. Since models including different extra-mixing pro- 
cesses predict distinct behaviours of Li vs. Be with age and 
with mass, the last few years have witnessed a flourishing of 
Be observations among field and cluster stars. The availability 
of state-of-the art high resolution spectrographs with high near- 
UV efficiency has made it possible to observe not only bright 
stars in close-by clusters, but also fainter members of more dis- 
tant, old clusters. 

These Be surveys have allowed constraining the mecha- 
nism responsible for the gap; more specifically, several stud- 
ies have now evidenced the presence of a correlation between 
Li and Be among cluster and field stars on the cool/red side of 
the gap (T e ff between ~ 6600 and 6300 K), providing strong 
support to rotationally induced mixing (e.g., Deliyannis et al. 
[TWgl Boesgaard et al. I2004al I2004bl l. As for cooler stars in 
older clusters, Randich et al. ( 120021 hereafter R02) found that 
solar-type members (T eff ~ 6150 - 5800 K) of M 67 and the 
2 Gyr old IC 465 1 have not depleted any Be and thus do not 
follow any Li-Be correlation. This finding allowed them to rule 
out a number of models; in particular, at variance with warmer 
stars on the red side of the gap, present models including rota- 
tionally induced mixing do not seem to well reproduce the ob- 
served Be vs. Li pattern for stars cooler than 6150 K in M 67: 
whereas they predict that at old ages these stars should have de- 
pleted some Be (a factor 2-3), no depletion is instead observed. 

In this work we present new Be measurements in three open 
clusters: the young NGC 2516 (150 Myr), the Hyades (600 
Myr), and the old M 67 (4.5 Gyr). In NGC 25 16 and the Hyades 
we mostly concentrated on cool MS stars (T e ff below 5400 K), 
the reason being that warmer stars in both the Hyades and other 
young clusters (such as the Pleiades and Alpha Per) have been 
extensively studied in the past by Boesgaard and collaborators 
(see Boesgaard et al. [2004b and references therein). In M 67 
we instead observed stars warmer than in the R02 sample, as 
well as two subgiants and two blue stragglers. Since the up- 
per turnoff of M 67 coincides with the cool edge of the Li gap, 
subgiants are supposed to be evolved from objects that were 



in the gap (Balachandran 1995); their Li and Be abundances 
are therefore important to constrain the gap evolution at old 
ages. As to blue stragglers, their light element content might 
provide insights into their formation mechanism. In summary, 
when joining our sample with available Be cluster data from 
the literature, we have a coverage of the T e ff interval from the 
Li gap down to 5100 K for three groups of stars, with ages of 
0.15, 0.6 and 4.5 Gyr. Based on this sample, our goal is to ob- 
tain a more complete empirical picture of the evolution of Be 
(and Be vs. Li) as a function of stellar age and temperature. 

2. Sample stars and observations 

The sample includes four members of NGC 2516, four mem- 
bers of the Hyades and eight members of M 67. The ob- 
servations were carried out in Service mode using UVES on 
VLT/Kueyen (Dekker et al. 2000) during the period November 
2001 - March 2002. All stars were observed using the same 
set-up; namely, UVES was operated in Dichroic Mode with 
Cross Dispersers #1 and #3 in the Blue and Red arms, respec- 
tively. Such a combination allowed us to cover the spectral 
ranges from ~ 3115 to 3940 A in the blue and from ~ 4780 
to 6810 A in the red. The 15fim pixels and the use of a 1 arcsec 
wide slit (projecting onto 4 pixels) and CCD binning 2x2 and 
1 x 1 in the blue and red yielded resolving powers R~ 40, 000 
and R~ 45, 000, respectively. The spectra were reduced by the 
ESO Quality Control Group using the UVES pipeline imple- 
mented in the MIDAS context. After careful inspection, we de- 
cided that the quality of the reduced spectra delivered by ESO 
did not require a new reduction of the raw frames; therefore, 
we just co-added spectra taken during different exposures, after 
applying the correction for radial velocity shifts and normal- 
ization. Total exposure times range between 3600 - 10800 s, 
depending on the star magnitude. The final S/N ratios per reso- 
lution element are in the intervals ~ 20 and 80 in the Be region 
and ~ 120 and 350 in the Li region. The sample stars together 
with the log of the observations and achieved S/N ratios are 
listed in TableQ] 

3. Abundance and Error Analysis 

Effective temperatures (T e ff's) for the majority of M 67 stars 
were retrieved from Jones et al. (1999 see also Randich et al. 
2006) who, in turn, had derived them from dereddened B-V 
colors using the calibration of Soderblom et al. (1993). We 
determined T e ff's for the two M 67 members not included in 
Jones et al. sample (S1034 and S1239) and for NGC 2516 
and Hyades stars in the same fashion. For M 67 a reddening 
E(B-V)=0.04 was used, following Jones et al., while we as- 
sumed E(B-V)=0.01 and 0.12 for the Hyades and NGC 2516. 
Adopted T e ff values are listed in the second column of Table[2] 
A surface gravity log g=4.5 was assumed for MS stars in 
NGC 2516 and the Hyades, while for M 67 we used the spec- 
troscopic values of Randich et al. (2006). Microturbulence 
was derived from log g and T e ft based on the relationship 
£ = 3.2x1 0- 4 (T e ff - 6390) - 1 .3(logg - 4. 1 6) + 1 .7 km/s (Nissen 
[T98TT) . 



S. Randich et al.: Beryllium in open clusters 
Table 1. Sample stars and log of the observations. Listed S/N ratios are per resolution element. 
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ID V (B-V)o Obs. Date Exp. time S/N @ Be S/N @ Li 

(sec) 



Main sequence stars 



M 67 S990 


13.43 


0.525 


Mar. 22,23 2002 


7200 


35 


150 


M 67 S995 


12.76 


0.521 


Feb. 3 2002 


5400 


40 


200 


M 67 S998 


13.06 


0.518 


Jan. 7 2002 


5400 


35 


175 


M 67 S2205 


13.14 


0.534 


Mar. 1,4 2002 


7200 


40 


200 


Hyades vB21 


9.14 


0.81 


Dec. 3 2001 


1200 


75 


350 


Hyades vB109 


9.40 


0.81 


Dec. 23 2001 


1800 


80 


350 


Hyades vB 182 


9.01 


0.83 


Dec. 21 2001 


1200 


80 


350 


Hyades vB187 


9.01 


0.75 


Nov. 29 2001 


1200 


80 


350 


NGC2516CTIO-2 


14.20 


0.79 


Dec. 19,23 2001 


10800 


20 


150 


NGC2516CTIO-3 


14.30 


0.82 


Dec. 20,22 2001 


10800 


20 


150 


NGC2516CTIO-6 


13.62 


0.66 


Dec. 30 2001 


9600 


30 


180 


NGC 2516DK213 


14.01 


0.76 


Dec. 28 2001 


5400 


20 


120 






Blue stragglers and subgiants 








M 67 S997 


12.13 


0.410 


Feb. 24 2002 


3600 


50 


210 


M 67 S2204 


12.89 


0.417 


Mar. 8, Apr. 6 2002 


5400 


40 


180 


M 67 S 1034 


12.65 


0.567 


Apr. 15,16 2002 


4800 


35 


200 


M67 S1239 


12.75 


0.692 


Mar. 7,8 2002 


5400 


30 


190 



Numbering for M 67 from Sanders ( 11977b 
Numbering for the Hyades from van Bueren Jl952l 
Numbering for NGC 2516 from Jeffries et al. dl998b 

Both Li and Be analyses were performed in the same way 
as in R02. Li abundances were derived by means of equivalent 
widths (EWs) of the Li i 6707.8 A feature and using Kurucz's 
WIDTH9 (Kurucz 1993) code and model atmospheres. Our 
spectral resolution, together with the low rotational velocity 
of all sample stars, allowed us to resolve the Li line from the 
nearby Fe i 6707.4 A feature and hence there was no need to 
correct measured Li EWs for its contribution. EWs and derived 
Li abundances are listed in Cols. 3 and 4 of Table[2] Errors asso- 
ciated to our Li abundances range between 0.08 and 0.17 dex. 
These combine the uncertainties in the EW measurement (as 
reported in Col. 3 of the table, and which include also the un- 
certainty due to the placement of the continuum) and in the 
stellar parameters (i.e. how much the Li abundance is affected 
by changes in the stellar parameters); as well known, the domi- 
nant source is the high sensitivity of Li to the effective temper- 
ature (0.08 dex per 100 K). Uncertainties in the stellar param- 
eters were taken from Randich et al (2006) for M 67 and are 
as follows: +70 K in T eff , +0.25 in log g, ±0.15 km/s in We 
assumed the same uncertainties for the Hyades and NGC 25 16. 
For a more detailed breakdown of all uncertainties affecting 
Li determination, we refer the reader to Charbonnel & Primas 
(120051 

As mentioned, Be abundances were determined using the 
same spectral synthesis recipes as in R02 because we would 
like to combine both samples for our final discussion. Namely, 
we synthesized a region of «5A, around the Be n doublet (A = 



3130.420 A and A = 3131.064 A), with Kurucz's model atmo- 
spheres and SYNTHE code (Kurucz 1993) and with the line 
list assembled and thoroughly tested by Primas et al. d!9971 l. 
Because of the severe blanketing present at these wavelengths, 
and because we are working with solar (or over-solar) metal- 
licity stars, the normalization of the spectra is critical. Since no 
line-free region is available around the Be n lines, we normal- 
ized our observed spectra by dividing them for the highest num- 
ber of counts in each spectrum, in a 15-20 A region centered on 
the Be n doublet; this step was performed before multiple spec- 
tra for a given object were co-added. The spectra were then 
further shifted (by a few %, on average) to match the computed 
synthetic spectra, but the overall continuum shape was main- 
tained. A few examples of spectrum syntheses are presented in 
Figs.Q~r|4] where we show our best-fit synthesis, two syntheses 
computed by varying the Be content by a fixed amount, and 
one synthesis computed without Be. 

Our final log n(Be) (as derived with ID model atmospheres 
and including a predicted Fe I line in our list of transitions) 
and associated uncertainties (see below) are listed in the last 
column of Table [2] 

As a final remark, we note that Randich et al. (2006} de- 
rived spectroscopic values of surface gravity for their sample 
stars in M 67 somewhat different from those assumed by R02; 
for this reason we carried out a new analysis of the two stars 
(S988 and S994) whose new log g values differed from the old 
ones by more than 0.25 dex (the assumed error in log g). New 
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Table 2. Effective temperatures, equivalent widths of the Li i 6707.8 A line, and derived Li and Be abundances. These are in the 
usual notation log n(X)=log N(X)/N(H)+12. Note that a positive offset of ~ 0.1 — 0.2 dex in log n(Be) might be present for stars 
cooler than 5400 K - see text. 



star 


T e ff 


EW(Li) 


log n(Li) 


log 


n(Be) 




(K) 


(mA) 










Main sequence stars 






M 67 S990 


6191 


45 ±4 


2.63+0.07 


0.96 


±0.13 


M 67 S995 


6210 


17 ±2 


2.10±0.08 


0.91 


±0.13 


M 67 S998 


6223 


44 ±2 


2.64±0.06 


0.91 


±0.13 


M 67 S2205 


6156 


11 ±3 


1.86+0.16 


0.81 


±0.13 


Hyades vB21 


5141 


< 3 


< 0.27 


0.85 


±0.21 


Hyades vB109 


5141 


5 ± 1 


0.49±0.11 


0.80 


±0.21 


Hyades vB182 


5079 


10 ± 1 


0.73±0.07 


0.80 


±0.21 


Hyades vB187 


5339 


17 ± 1 


1.27+0.06 


0.80 


±0.21 


NGC2516CTI0-2 


5238 


147 ± 3 


2.69+0.06 


1.16 


±0.21 


NGC 2516CTI0-3 


5110 


174 ±4 


2.80±0.07 


1.04 


±0.21 


NGC2516CTIO-6 


5659 


124 ±5 


2.93±0.07 


1.21 


±0.21 


NGC 2516DK213 


5305 


108 ±2 


2.41+0.06 


1.06 


±0.21 




Blue stragj 


*lers and subgiants 






M 67 S997 


6700 


< 5 


< 1.91 


0.31 


±0.25 


M 67 S2204 


6668 


<5 


< 1.89 


-0.20 


±0.29 


M67 S1034 


5969 


< 3 


< 1.14 


0.12 


±0.13 


M 67 S 1239 


5477 


8± 1 


1.12+0.10 


0.11 


±0.21 



and old abundances for these stars are listed in Table [3] In the 
following we will consider the new abundances derived in this 
work. 

3.1. Errors in Be abundances 
3.1.1. Random uncertainties 

In order to estimate the uncertainty associated to our Be abun- 
dances, we computed the sensitivity of the derived Be abun- 
dances to changes of the stellar parameters and of the place- 
ment of the continuum. Our final errors do not include the un- 
certainty related to the atomic physics of the Be transitions (the 
log gf values are supposedly known with high accuracy) or the 
absolute error on the model atmospheres, which is difficult to 
properly estimate. Because our sample includes stars in dif- 
ferent evolutionary stages, we performed our sensitivity analy- 
sis for three stars characterized by different stellar parameters: 
S998, vB182, and S2204. We chose S998 as representative of 
our MS M 67 sample, vB182 as representative of the cooler 
MS stars, and S2204 as one of the hottest stars in our sample 
with very little Be left. Because the metallicity of these clus- 
ters is rather well determined, we assumed an uncertainty of 
+0.05 dex on [Fe/H]. For the placement of the continuum, we 
estimated that +3% is representative of this uncertainty. This 
number was derived by normalizing our observed spectra with 
a different method, i.e. using the computed synthetic spectra as 
the main reference to be matched. A summary of our Be error 



Table 3. New (present analysis) and old (R02) beryllium abun- 
dances for S988 and S994 already analyzed by R02. 



star T eff log g oW log n(Be) | d log g new log n(Be), 

S988 6153 4.5 0.88 4.1 0.65 

S994 6151 4.5 1.16 4.1 0.86 



analysis is provided in Table [H where the numbers reported in 
the last row (<x tot ) are the square root of the sum in quadrature 
of all factors. We note that when the dependence on some stel- 
lar parameters was not perfectly symmetric (when increasing 
or decreasing that parameter by the chosen fixed amount), we 
always took the average as the representative value. 

3.1.2. Systematic uncertainties 

We remind from R02 that the list of lines we have used includes 
a predicted line (i.e., not fully identified in the laboratory) of 
neutral iron at A — 3131.043 A which is needed in order to ob- 
tain a satisfactory fit of the A = 3131.064 A Be line in the solar 
spectrum. Our analysis of the Sun yields a photospheric solar 
Be abundance log n(Be) = 1.11, in excellent agreement with 
the value found by Chmielewski et al. d!9751 l. The atomic char- 
acteristics of this line were constrained by Primas et al. ([1997 ) 
on a large sample of stars, spanning a wide range of stellar pa- 
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Table 4. Sensitivity of Be abundances to T eff , log g, [Fe/H], and placement of the continuum, for three stars of our sample. 



S998 S2204 vB182 



Teff(K) 


6222 


6668 


5079 


±70 K 


±0.03 dex 


±0.07 dex 


±0.05 dex 


logg 


4.0 


4.0 


4.5 


±0.25 


±0.10 dex 


±0.08 dex 


±0.17 dex 


f (kms- 1 ) 


1.50 


1.35 


0.84 


±0.15 km s" 1 


0.00 


0.00 


±0.05 dex 


[Fe/H] (dex) 


0.00 


0.00 


0.13 


±0.05 dex 


+0.03 dex 


+0.10 dex 


+0.05 dex 


±3% cont. 


+0.04 dex 


+0.25 dex 


+0.08 dex 


0" tot 


±0.13 dex 


±0.29 dex 


±0.21 dex 



rameters (in temperature, gravity and metallicity). It is impor- 
tant to note, however, that Garcia Lopez et al. (|T995) proposed 
an alternative solution to the unsatisfactory fit of the redder 
component of the Be doublet, i.e. an increase of the oscillator 
strength of the Mn I line at X = 3131.037 A by +1.5 dex (from 
the value of - 1 .725 reported in the Kurucz linelist). They found 
that such line becomes the dominant blending transition of the 
redder component of the Be doublet for stars with effective 
temperature below 5200 K, which prevented them from making 
a full detection of Be in their lowest temperature stars. Because 
of this finding and of the fact that, at solar (or above solar) 
metallicity, the predicted line used by Primas et al. ( 1997) was 
tested only down to ~ 5400 K (and our current sample reaches 
instead 5079 K), we decided to run a few more tests and to 
compare the effect of the two proposed solutions, especially 
for the lowest temperature stars. 



First of all, we confirmed that both solutions do not have 
any important effect for stars above 5400 K, except for allow- 
ing a better fit of the blue wing of the A = 3131.066 A. As we 
tested lower effective temperatures, we found that the two so- 
lutions differ by 0.2 dex in the final derived Be abundance in 
the 5000-5100 K range, while in the 5250-5400 K interval the 
difference amounts to 0.15-0.10 dex, with our solution always 
yielding lower Be abundances. In other words, for stars below 
5400 K the strength of the predicted Fe I line is larger than that 
of the the Mn I and the difference increases with decreasing 
temperature. Note, however, that our proposed solution still al- 
lows us to measure Be abundances down to the lowest T e jf in 
our sample. Whereas there is no solid argument in favour of one 
solution with respect to the other, depending on the solution 
one decides to follow, the lowest temperature stars may end up 
having slightly different abundances, with the offset depending 
on the temperature. For this reason we are not able to derive a 
definitive absolute Be abundance and thus to quantify the ex- 
act amount of Be depletion suffered by Hyades members in our 
sample; nevertheless the relative abundances of NGC 2516 and 
Hyades stars with similar temperatures suggest that the latter 
have undergone some depletion. 



4. Results 

4.1. Comparison with previous studies 

Our Li abundances for most stars are in very good agreement 
with the values derived in previous studies (Jeffries et al. 1998 
for NGC 2516; Thorburn et al. [T9931 for the Hyades; Jones 
et al. 119991 and Balachandran 119951 for M 67). For one star 
only, CTIO-6 in NGC 2516, we derive an abundance some- 
what lower than (but given the errors still consistent with) the 
value of Jeffries et al. (1998); namely, we obtain log n(Li)=2.93 
to be compared with value log n(Li)=3.04. Note however that 
Jeffries et al. (1998) measured a much higher EW (212 mA), 
which in our abundance scale would correspond to log n(Li) ~ 
3.4. Also, we detected the Li line in the spectra of two stars 
for which only upper limits were available in the literature 
(S2205: log n(Li) pres =1.86 vs. log n(Li)u t < 2.34 -Jones et 
al. [19991 S12 39: lo g n(Li) pres =1.12 vs. log n(Li) lit . < 1.1 - 
Balachandran 1995) and were able to put more stringent upper 
limits on Li abundances for three stars: namely S997 and S2204 
(log n(Li) pres . < 1.91 vs. log n(Li) Ht . < 2.46; log n(Li) pres . < 
1.89 vs. log n(Li) Ht < 2.36 -Shetrone & Sandquist [20001 
and S1034 (log n(Li) pres . < 1.14 vs. log n(Li)i it . < 1.6 - 
Balachandran [T9"9"5l 

Our study provides the first determination of Be for 
all the sample stars, but vB21; this star was included in 
the study of Garcia Lopez et al. i 19951 1 who measured 
log n(Be)=0. 90+0.25, in agreement with our own determina- 
tion (log n(Be)= 0.85 + 0.21). 

4.2. Be vs. T e ff 

In Fig. [5] we plot log n(Be) (upper panel) and log n(Li) (lower 
panel) as a function of T e fr for our sample stars and the sam- 
ple of R02. The present M 67 sample includes four new stars 
warmer than 6150 K, but no cooler cluster members; therefore 
the conclusions of R02 for solar-type stars remain unchanged: 
independently of their absolute Be abundance, almost all the 
stars in the range 5600 < T e ff < 6150 K show, within the errors, 
the same relative abundance, in spite of the different levels of Li 
depletion (see lower panel of Fig. |5j. On the other hand, stars 
warmer than 6150 K both in the present and in the R02 sam- 
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pies have lower abundances, suggesting that they have suffered 
a certain amount of Be depletion. Note that, as discussed in de- 
tail in Sect. 15. II below, we consider the abundance of SHJM2, 
a member of the 55 Myr old cluster IC 2391, as representative 
of the value of the initial Be abundance in our scale. 

As to stars cooler than 5400 K -note that our sample does 
not include stars with T e fr in the range 5400-5600 K-, Be abun- 
dances of NGC 25 16 members are consistent with no depletion 
in spite of the factor ~2-8 Li depletion. The significantly more 
Li-poor Hyades have lower Be abundances: the average value 
of the four Hyades stars is ~ 0.25 dex (or a factor 1.8) below 
the average of the three NGC 2516 stars cooler than 5400 K. 
This result would not change if the abundances were derived 
assuming the solution proposed by Garcia Lopez et al. (119951 
-see discussion in Sect. [3]). 

Finally, both the two blue stragglers (S997 and S2204) and 
the two evolved M 67 members are characterized by significant 
amounts of Li and Be depletion. 

5. Discussion 

5.1. Initial beryllium abundance 

As discussed by R02, the Be abundance of the young star 
SHJM2 in IC 2391 (log n(Be)=l.ll), as well as our value for 
the photospheric solar Be (log n(Be)=l.ll), are about a factor 
of two lower than the meteoritic abundance (log n(Be) = 1 .42 - 
Anders & Grevesse 1989LU). Abundances for stars in the present 



study and in R02 are all below this value. Also, maximum Be 
abundances similar to ours are found by Garcia Lopez et al. 
d 19951 1 for the Hyades and by Santos et al. (2004) for their so- 
lar metallicity stars; both studies derive a Be abundance for the 
Sun very close to our determination. IC 2391 has a solar metal- 
licity (Randich et al. 120011) and thus its initial Be abundance 
is most likely the same as for the solar system. The difference 
between our log n(Be) values for the Sun and SHJM2 and the 
meteoritic Be abundance suggests that either a) both the Sun 
and SHJM2 have depleted a factor of two Be; or b) they have 
not undergone any depletion, the difference with the meteoritic 
value is due to systematic effects in the adopted analysis, and 
the log n(Be) value of SHJM2 is indicative of the initial Be 
abundance in our scale. We mention in passing that the discrep- 
ancy between the solar photospheric Be and meteoritic abun- 
dance has been known since long time and it has been discussed 
in several works. 

Hypothesis a) and, specifically, the fact that SHJM2 has 
depleted some Be, would imply that Be burning occurs dur- 
ing the pre-main sequence (PMS) phases, at variance with the 
observational evidence of no PMS Li depletion for stars in this 
T e jf range (Randich et al. 2001 and references therein); it would 
indeed be hard to explain how this star depleted some Be while 
maintaining its initial Li. For this reason, we regard possibil- 
ity a) as very unlikely and favor hypothesis b) above. At the 
same time, since the metallicities of M 67 and NGC 2516 are 
also solar (Randich et al. 120061 Jeffries et al. 1998 ), it is un- 
likely that these clusters had different initial Be abundances due 



1 A slightly lower value, log n(Be) = 1.38 + 0.08, has been recently 
reported by Asplund et al. ( 2005 ). 



to Galactic evolutionary effects. We conclude therefore that all 
the sample stars with log n(Be) consistent with the abundance 
of SHJM2, including the Sun, have most likely not undergone 
any Be depletion. A strong support to this hypothesis is given 
by the analysis of Balachandran & Bell ( 1998 ) of the solar Be 
who showed that, taking into account an additional amount of 
UV opacity, needed to solve the inconsistency between solar 
oxygen abundance derived from IR and UV lines, accounts per- 
fectly for the missing Be and makes the discrepancy between 
the solar photospheric Be and the meteoritic value disappear. 
The conclusion that the photospheric solar Be is consistent with 
the meteoritic one was also reached by Asplund et al. (2004); 
our assumption hence appears fully justified. We finally men- 
tion that a higher value of the solar Be abundances was derived 
by Boesgaard et al. d2003l log n(Be) G = 1.30), by using 2002 
version of MOOG which includes all of the Kurucz opacity 
edge; also, for undepleted cluster stars in their sample, they de- 
rive maximum abundances higher than ours (and in particular 
than that of SHJM2). These differences suggest the presence of 
an offset between the abundance scales. 



5.2. Be and Li depletion 

In Fig. [6] we plot log n(Be) as a function of log n(Li) for 
our sample stars, the sample of R02, and the Hyades (from 
Boesgaard et al. 2004a and Garcia Lopez et al. 1995 ). In the up- 
per, middle, and lower panels we show stars in three tempera- 
ture ranges: namely, 6300>T eff > 6150 K, 5600 <T eff < 6150 K, 
and T e fr < 5400 K. This subdivision broadly corresponds to the 
different regimes of Be depletion seen in Fig. [5J specifically, 
stars that show both Li and Be depletion, stars that show Li 
depletion, but no Be depletion, stars that show Li depletion 
and may show Be depletion. Whereas the precise boundaries 
between these three intervals might be slightly different, it is 
important to subdivide the sample stars in subgroups since, as 
we discuss below, they likely correspond to different dominant 
mixing processes. In the upper panel we show for comparison 
also Hyades stars on the red side of the Li gap. We finally note 
that our sample does not include stars with T e ff between 5400 
and 5600 ; this T e ff range is thus not represented in the figure. 

The figure evidences distinct behaviors for stars in the three 
subsamples. Namely, the middle panel, which includes only 
stars from R02 plus the Hyades, clearly indicates the absence of 
a Be vs. Li correlation for solar-type stars in the 6150-5600 K 
range in both old and young clusters: they span two orders 
of magnitude in Li abundances, they have different ages and 
metallicities, but they share the same Be content. This in turn 
implies that the mixing mechanism responsible for MS Li de- 
pletion in this temperature range does not extend deep enough 
to cause also Be depletion. On the other hand, the enlarged 
sample with respect to R02 shows that a correlation between Li 
and Be abundances is present for warmer stars in M 67; these 
stars rather well fit into the Be vs. Li trend for the Hyades on 
the red side of the Li gap. 

Finally, the bottom panel and, in particular the difference 
between Be abundances of the Hyades and NGC 2516 mem- 
bers, suggest that in the Hyades cool stars Li depletion is ac- 
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companied by some amount of Be depletion. It would be not 
formally correct to claim a Be-Li correlation for these stars, 
since NGC 25 16 members cover a large range in Li abundance, 
but have not depleted any Be; also, Hyades stars themselves 
show different amounts of Li depletion, but the same Be con- 
tent. Nevertheless, the figure shows that, in this T e ff regime, at 
variance with solar-type stars, but similarly to warmer objects 
(T e ff > 6 1 50 K), the mechanism responsible for Li depletion can 
reach deep enough to cause also Be depletion. We stress again 
that, whereas the absolute Be abundances of the coolest sample 
stars in the figure would be higher had we used the solution of 
Garcia Lopez et al. d 19951 1. the relative depletion would remain 
the same. 

5.3. Cool stars: Be as a function of age 

Garcia Lopez et al. ( 1 19951 1 already found evidence of Be deple- 
tion for two cool Hyades members, while Santos et al. (f2004) 
detected Be depletion in several old field stars cooler than about 
5600 K. Our results not only reinforce the conclusion that cool 
stars might deplete Be, but most importantly, the comparison 
between stars with similar temperatures and different ages al- 
lows us to put constraints on the age when Be depletion starts 
on the MS. 

In Fig. [7] we plot log n(Be) as a function of T e ff for 
NGC 2516, the Hyades (present, Garcia Lopez et al. [1995l and 
Boesgaard et al. 2004a samples), the Pleiades from Boesgaard 
et al. (120031) . and old field stars from Santos et al. (12004b . Note 
that the original sample of Santos et al. is mostly composed 
by stars older than 1 Gyr and we do not plot in the figure the 
four young stars in their sample. The figure shows that Be de- 
pletion starts being present below ~ 5600 K: for T e fr between 
5400 and 5600 K Be depletion is seen only among old field 
stars, while for temperatures below 5400 K, Be depletion is 
seen also among the Hyades, but not among NGC 2516 mem- 
bers. This implies that the age at which Be depletion begins in 
cool stars increases with increasing T e ff: namely, stars cooler 
than ~ 5400 K start depleting Be between 150 and 600 Myr, 
depletion starts after 600 Myr for stars in the 5400-5600 inter- 
val, while stars warmer than 5600 K (but cooler than 6150 K), 
including the Sun, never deplete Be. 

Li depletion in MS stars cooler than 5500-5400 K is nor- 
mally explained with convection, since the convective enve- 
lope reaches deep enough to mix Li-poor material into the sur- 
face. The question then arises whether Be depletion is also ex- 
pected for these stars based on standard models. In order to 
address this question, we calculated evolutionary tracks for a 
0.9 M & star, corresponding to the mass of our coolest stars 
in the Hyades (see below), using the models by Sestito et al. 
(2006). These were computed employing an updated version of 
the evolutionary code FRANEC (e.g. Chieffi & Straniero [T9"89l > 
and adopting updated physical inputs (equation of state, opac- 
ity tables, nuclear cross sections, etc., see Sestito et al. 2006 
for further details). As customary in standard models, the only 
mixing mechanism included is convection, whose extension 
was determined by adopting the classical Schwarzschild cri- 
terion and the mixing length formalism. Element diffusion has 



been included during the MS (Ciacio et al. 1997 with diffusion 
coefficients from Thoul et al. II 9941) . Three different models 
were computed; namely, we adopted two different solar mix- 
tures: the solar composition by Asplund et al. (2004; hereafter 
As04) and that by Grevesse & Noels ( [19931 hereafter GN93). 
The first track computed is the "standard" model by Sestito et 
al. ( |2006l model 1 in their Table 1): As04 mixture, Z = 0.013 
(corresponding to [Fe/H]=0 for this composition), He content 
Y = 0.27, and mixing length parameter a = 1.9. The other two 
models have the same He abundance and a but the Fe content of 
the Hyades ([Fe/H]= +0.13) has been considered, in one case 
adopting the As04 mixture (and the resulting global metallicity 
is Z = 0.018), while in the other using the GN93 composition 
(Z = 0.025). We mention that at the age of the Hyades, the 
[As04, Z = 0.013] model has an effective temperature of 5400 
K, while the [As04, Z = 0.018] model has T eff =5200 K and the 
[GN93, Z = 0.025] one has T eff =5130 K, consistent with the 
temperatures of Hyades stars in our sample. 

In Fig. [8] we plot the temperature at the base of the convec- 
tive envelope (Tcz) predicted by our models as a function of 
time: the figure clearly shows that, for all the models, Tcz is 
slightly higher than the Be ignition temperature only for ages 
younger than about 20 Myr, implying that Be depletion can in 
principle occur only during PMS, while no Be depletion is ex- 
pected at older ages. This is shown in Fig. [9] where we plot the 
resulting evolution of Be/Beo as a function of age: consistently 
with the discussion above, none of the three models predict 
significant Be depletion. The model with solar metallicity pre- 
dicts a slightly lower amount of Be depletion with respect to 
the other two, as expected (the higher is Z, the higher is the 
efficiency of convection). However, and most importantly, the 
difference between the solar metallicity track and the two mod- 
els with the Fe content of the Hyades is almost negligible, and 
in all cases a very small amount of Be is depleted with respect 
to the initial value. In particular, we note that even PMS Be 
depletion is always very small (almost negligible), since Tcz 
reaches values only slightly higher than the temperature of Be 
ignition and, furthermore, the exact value of ignition tempera- 
ture depends on the physical assumptions in the code. We men- 
tion in passing that this result is in agreement with the early 
calculations of Bodenheimer ( 1966 ) 

Our conclusion is that Be depletion observed in the cool 
Hyades with respect to NGC 2516 cannot be explained only 
by convective mixing and cannot be reproduced in the frame- 
work of standard models. An additional, or extra-mixing mech- 
anism must be present, analogously to what invoked to explain 
Li depletion in more massive, warmer stars. Most obviously, 
the mechanism is not necessarily the same. Our Be determi- 
nations suggest that this physical process starts being efficient 
after an age of ~ 150 Myr. 

5.4. Stars out of the main sequence 
5.4.1 . Subgiants 

M 67 subgiants SI 034 and S1239 studied here were also in- 
cluded in the Li survey of Balachandran (1995), who de- 
rived for the two stars masses and ZAMS temperatures equal 
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to 1.282 M and 6349 K (S1034) and 1.3 M and 6450 K 
(S1239); as argued by Balachandran, their low Li abundances 
strongly support the hypothesis that they have depleted Li 
during the MS phases and have evolved from the Li gap re- 
gion. Our study provides the first measurements of Be in stars 
evolved from the gap, thus allowing us to carry out further com- 
parison with model predictions. Both stars show similar factors 
of Be depletion (a factor of about 10 -see Table |2), suggest- 
ing the presence of a Be gap, similarly to what found for the 
Hyades MS stars Boesgaard & King (2002). As discussed at 
length in the series of papers by Boesgaard and collaborators 
and mentioned in Sect. Q] the tight correlation between Be and 
Li depletion for stars in the gap suggests that light element de- 
pletion is driven by rotational mixing. Our Be determinations 
for the two subgiants are consistent with this hypothesis. 

More quantitatively, Sills & Deliyannis (2000) presented 
evolutionary models of subgiants in M 67 including three mix- 
ing processes (diffusion, mass loss, rotational mixing) and fol- 
lowing the evolution of both Li and Be abundances. They con- 
cluded that Li measurements are in good agreement with mod- 
els including mixing driven by rotation; they stressed however 
that Be data for subgiants in M 67 would also be very helpful 
in order to constrain the various free parameters entering the 
models. Our measurements of Be for two subgiants in M 67 
represent a first step towards that direction. In Fig.[l0]we com- 
pare our results with the predictions of the model of Sills & 
Deliyannis including rotationally induced mixing; while the 
other models discussed in that paper are not shown, since they 
were already ruled out the basis of Li alone. Note that the 
T e ff sequence in the figure represents an evolutionary sequence 
for a given mass (1.26 M -somewhat below the value given 
by Balachandran 1995 for subgiants in M 67) rather than a se- 
quence of different masses. The figure evidences a very good 
agreement between the model and our measurement for S 1239, 
the cooler, slightly more massive and more evolved star, whose 
Be content is the result of MS depletion plus post-MS dilu- 
tion. The agreement is considerably worse for the less evolved 
SI 034, whose Be content is the result of MS depletion only, 
and which appears more Be depleted than predicted (but also 
its Li abundance is below model predictions -see Fig. 4 in 
Sills & Deliyannis). We suggest that the reason for this dis- 
crepancy could be an initial higher rotational velocity for this 
star than the one assumed by Sills & Deliyannis in their model 
(Vm. = 10 km/s), which would result in a larger amount of light 
element depletion on the MS. 



5.4.2. Blue stragglers 

Light elements also provide an useful tool to investigate the 
mechanism of formation of blue stragglers and, in particular, 
whether they are the product of collision processes or binary 
mergers. To our knowledge, no quantitative models exist to- 
date presenting predictions for Li and Be abundances for the 
different scenarii of blue straggler formation. However, Sills et 
al. ([T997 ) qualitatively suggested that in the case a blue strag- 
gler is the result of a binary merger, it should not have observ- 



able Li; on the other hand, collisions would not result in any Li 
(and a fortiori Be) destruction. 

The first observations of light elements in blue strag- 
glers belonging to M 67 were performed by Pritchet and 
Glaspey (1991), who measured Li in seven stars; they did not 
detect the Li line in any of them and concluded that some form 
of mixing, possibly due to binary coalescence, had affected the 
outer envelope of blue stragglers. We note however, that up- 
per limits for most of their sample stars were indeed consistent 
with abundances of Li-poor MS cluster stars. 

In a more recent study, Shetrone & Sandquist (2000| car- 
ried out new Li measurements of blue stragglers in M 67, in- 
cluding both S997 and S2204, the two stars studied here. We 
refer the reader to Shetrone & Sandquist (2000) for a detailed 
discussion of the properties of the two stars; we just recall that 
the high eccentricity binary nature of one of the two (S997) 
suggests that it is unlikely a binary merger, which is instead the 
most plausible scenario for the formation of S2204. Shetrone 
& Sandquist (2000) did not detect the Li line in either star; 
although the non-detection of Li might favour the merger sce- 
nario, these authors point out that their upper limits eventu- 
ally do not provide any clues on the formation mechanism. The 
upper limits are not very stringent and are consistent with the 
abundances of MS unevolved cluster stars (i.e., the presence 
of Li cannot be ruled out); furthermore, even a low Li abun- 
dance could be explained by the fact that the two stars fall in 
the Li gap T e jj interval and thus could have undergone signifi- 
cant depletion after the blue straggler is formed in the collision 
scenario. 

Until the present study no Be determinations for blue strag- 
glers in open clusters were available in the literature. We have 
now determined Be abundance for S997 and S2204, finding 
that they are both depleted, but affected by different amounts 
of Be depletion; namely, a factor of about 6 and 20 for S997 
and S2204, respectively. At variance with Li, the low Be abun- 
dances are not consistent with those of MS stars in M 67; as 
in the case of Li, we cannot completely exclude that the two 
stars have depleted Be after formation. Nevertheless, we no- 
tice that Be abundance for star S2204 is much below the values 
measured for stars in the Be gap. Thus, whereas Be abundance 
in S997 would not be inconsistent with either scenario, we sug- 
gest that the factor of three higher depletion measured in S2204 
might be an indication of a different formation of the two sys- 
tems, with this star being the product of a binary merger. 

6. Conclusions 

We have carried out Be observations in three open clusters of 
different ages: NGC 2516, the Hyades, and M 67. The present 
study, together with available data from R02 for M 67 and from 
other previous Be surveys of open clusters, allows us to cover 
a grid of MS stars which span from 150 Myr to 4.5 Gyr in age 
and from the Li gap down to 5100 K in effective temperature. 

We confirm that stars in different T e g interval are character- 
ized by different amounts of Be depletion and Be vs. Li behav- 
iors; namely, i. stars warmer than 6150 K are characterized by 
some Be depletion and show a Be vs. Li correlation; ii. Li and 
Be are instead not correlated for solar-type stars. Whereas Li 
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depletion ranges between null (for the young star in IC 2391) 
and a factor larger than 100 (in the Sun and the Li poor stars 
in M 67), none of the stars show significant Be depletion; iii. 
the Hyades stars cooler than 5400 K are a factor of ~ 1.8 more 
Li depleted than stars of similar temperature in the younger 
NGC 2516; iv. subgiants and blue stragglers are severely de- 
pleted in Be. 

The result for warm stars is consistent with the extensive 
work carried out by Boesgaard and collaborators and, in par- 
ticular, with the observed Li-Be correlation for stars in the Li 
gap in the Hyades and other clusters of similar age. As men- 
tioned in Sect. [1] the presence and slope of this correlation are 
in agreement with the predictions of models including slow 
mixing driven by rotation. Be depletion in the two subgiants 
in M 67 that are evolved from stars in the gap is also consis- 
tent with rotationally induced mixing and our measurements 
actually can provide additional constraints on theoretical mod- 
els. On the contrary, the lack of a correlation between Li and 
Be for solar-type stars implies that the main process at work in 
these stars must be different. R02 indeed ruled out models in- 
cluding mixing induced by rotation and diffusion which predict 
a Be-Li correlation. They noticed that models including grav- 
ity waves would well fit the observed Be vs. Li distribution, but 
those models would not be able to explain the observed scatter 
in Li in M 67. 

In a more realistic approach one might imagine that var- 
ious processes (rotationally driven mixing, meridional circu- 
lation, diffusion, gravity waves, tachocline) are all present at 
the same time; different mechanisms might then either become 
dominant or counter-balance in different temperature regimes. 
In this context, very recently Charbonnel & Talon (2005) pre- 
sented a model including both mixing driven by rotation and 
gravity waves. They showed that this model is able to repro- 
duce the solar Li abundance as well as the internal rotational 
profile of the Sun and that models with different initial rota- 
tional velocities would result in a range of Li abundances, con- 
sistent with the observed scatter in M 67. Since Charbonnel & 
Talon (120051 ) did not provide any predictions for the evolution 
of Be abundances, we cannot make any quantitative compari- 
son with the results of the present study. We also mention in 
passing that their model would not be able to reproduce the 
plateau in Li abundances that is observed for ages greater than 
~ 2 Gyr (Randich 2006 and references therein). More work is 
therefore needed on theoretical grounds; we stress however that 
considering the interplay between different processes might in- 
deed represent the way towards the solution to the light element 
depletion problem. 

As to cooler stars, whereas between 5600 K and 5400 K Be 
depletion is seen only in field stars older than 1 Gyr, when we 
reach stars as cool as 5400 K, Be depletion starts at younger 
ages, between 150 and 600 Myr. We have shown that this can- 
not be explained by the effect of increasing depth of the convec- 
tion zone (that instead explains Li depletion), and that an addi- 
tional source of extra-mixing is needed also for cool stars. No 
models including extra-mixing processes are so far available 
for stars significantly cooler than the Sun; our results clearly 
point toward the necessity of developing such models. 



Finally, our Be determination for the two M 67 blue strag- 
glers does not allow us to draw any definitive conclusion on 
their formation scenario; nevertheless, the significant differ- 
ence in their Be content suggests that they have formed by two 
different processes. 
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Fig. 1. The observed spectrum of the M 67 star S998 in the Be region with four spectrum syntheses overplotted: our best-fit 
synthesis (continuous line), two syntheses computed respectively with a variation of the Be content of +0.15 dex (dashed lines) 
and one synthesis computed without beryllium (dotted line). The two vertical bars identify the components of the Be n doublet. 
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Fig. 2. Same as Fig.Q] but for the M 67 star S2204 and for a variation in the Be content of +0.3 dex. 
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Fig. 3. Same as Fig.Q] but for the Hyades star vB182 and for a variation in the Be content of ±0.25 dex. 
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Fig. 4. Same as Fig.Q] but for the NGC 2516 star DK213 and for a variation in the Be content of +0.25 dex. 
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Fig. 5. Upper panel: Beryllium abundances (log n(Be), in the usual logarithmic scale where log n(H)=12) as a function effective 
temperature (T e ff) for our sample stars plus the sample of R02. Symbols are as follows: M 67: filled circles (present sample) and 
crossed filled circles (R02); NGC 2516: open squares; Hyades: open circles; IC 4651 (R02): open diamonds; 1C 2391 (R02): 
open star symbol. Lower panel: same as upper panel, but Li -log n(Li)- abundances are plotted. Names of M 67 stars out of the 
MS are indicated in both panels. Horizontal lines denote initial Be and Li abundances in our scales. 
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Fig. 6. log n(Be) as a function of log n(Li) for our sample stars, the sample of R02, and the Hy ades from Boesgaard et al. ( 2004a ) 
and Garcia Lopez et al. ([1995 ). Upper panel: stars warmer than 6150 K. Middle panel: stars with 5600 < T e a < 6150 K; Lower 
panel: stars cooler than 5400 K. M 67 members out of the MS are not plotted in the figure. The Sun is also shown in the middle 
panel. Symbols are the same as in Fig.[5j Hyades from the literature, not plotted in Fig [5] are represented here as open circles 
with stars (in all panels for T e g < 6300 K) and open triangles (upper panel, stars with T e g 6300 K). Horizontal and vertical dashed 
lines denote initial Be and Li abundances in our scales. 
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Fig. 7. log n(Be) as a function of T e ff for NGC 2516 (150 Myr -open squares), the Hyades (600 Myr -present sample: open 
circles; sample of Boesgaard et al. 2004a open circles with stars; sample of Garcia Lopez et al. 119951 open circles with plus 
symbols), Pleiades (125 Myr -open triangles; data of Boesgaard et al. 120031 . and field stars (> 1 Gyr -filled triangles; sample of 
Santos et al. 120041 . The horizontal line denotes the initial Be abundance in our scale. 
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Fig. 8. Model predictions for the temperature at the base of the convective zone as a function of age. The solid horizontal line 
indicates the temperature necessary for Be to burn; the solid, dashed, and dotted curves show the predictions of the three different 
models described in the text. 
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Fig. 9. Model predictions for the evolution of N(Be)/N(Be)o as a function of age for the three different models. 
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Fig. 10. Comparison of model predictions of log n(Be) for subgiant stars in M 67 with our measurements for S1034 and S1239. 
The figure is based on Figure 6 in Sills & Deliyannis (2000), but the initial Be abundance has been scaled to our value. With 
respect to that figure we show only the model including rotationally driven mixing. 



